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A new algorithm for calculating the spin- and parity-dependent shell model nuclear level densities 
using the moments method in the proton-neutron formalism is presented. A new, parallelized code 
based on this algorithm was developed and tested using up to 4000 cores for a set of nuclei from the 
sd-, pf-, and pf + gg/2 - model spaces. By comparing the nuclear level densities at low excitation 
energy for a given nucleus calculated in two model spaces, such as pf and pf + gg/2, one can 
estimate the ground state energy in the larger model space, which is not accessible to direct shell 
model calculations due to the unmanageable dimension. Examples for the ground state energies of 
for ^''Ge and ^^Se in the pf + gg/2 model space are presented. 

PACS numbers: 21.10.Ma, 21.10.Dr, 21.10.Hw, 21.60.Cs 



I. INTRODUCTION 

Spin- and parity-dependent nuclear level densities rep- 
resent an important ingredient for the theory of nuclear 
reactions. For example, the Hauser-Feschbach approach 
[H requires exact knowledge of nuclear level densities for 
certain quantum numbers J'^ of spin and parity in the 
Gamow window of excitation energies around the parti- 
cle threshold In most of the cases relevant to nu- 
clear astrophysics, where experimental information is not 
available, the reaction rates for medium and heavy nu- 
clei can only be estimated using the Hauser-Feshbach ap- 
proach. Nuclear level densities are usually obtained using 
the back-shifted Fermi gas approximation 0-01, which 
was improved over the years. More modern approaches 
to the level densities based on the mean field were re- 
cently proposed by Goriely and collaborators 0, S] , and 
0. These approximations assume an independent parti- 
cle model in a mean field that lacks information about the 
many-body correlations. These correlations can be in- 
cluded exactly if one can fully diagonalize the many-body 
nuclear Hamiltonian, a task of increasing difRcultj^ Al- 
ternatively, one can use Monte- Carlo techniques 
or other methods of the statistical spectroscopy 
including applications to large shell-model spaces 

Most of these methods 0, [13, [H calculate the 

density of states and later use a spin-weight factor that 
includes an energy-dependent cut-off parameter to ex- 
tract the spin-dependent nuclear level density. Although 
there are recent efforts to improve the accuracy of such 
parametrizations [2l| . It was shown that the spin cut-off 
parameter has very large fluctuations at low excitation 
energy, when compared with the shell model results [12] ■ 
Statistical spectroscopy provides a path to a direct cal- 
culation of the spin cut-off parameter using a polynomial 
expansion for its estimate 23] . This approach was re- 
cently investigated (see Fig. 2 of Ref. [13]), where it was 
shown that although the smooth part of the energy de- 
pendence of the spin cut-off parameter can be described 
reasonably well, significant fiuctuations are still present 




in the low-energy regime. The quality of the results of 
this approach are mixed. Therefore, one would like to 
have a spin-projected method of calculating nuclear level 
densities that is accurate and fast. The parity is usu- 
ally taken as equally distributed, although there are at- 
tempts [3, [13] to model the effect of the uneven parity- 
dependence of the level densities at the low excitation 
energies of interest for nuclear astrophysics. 



Recently, we developed a strategy [22, 125l - l27| of calcu- 
lating the spin- and parity-dependent shell model level 
density. The main ingredients are: (i) extension of meth- 
ods of statistical spectroscopy [H, [13] by exactly calcu- 
lating the first and second moments for different configu- 
rations at fixed spin and parity; (ii) exact decomposition 
of many-body configurational space into classes corre- 
sponding to different parities and number of harmonic 
oscillator excitations; (iii) development of new effective 
interactions for model spaces of interest starting with the 
G-matrix [29| and fixing/fitting monopole terms and/or 
linear combinations of two-body matrix elements to ex- 
perimental data; and (iv) an accurate estimate of the 
shell model ground state (g.s.) energy. The calculation 
of the latest ingredient is generally as time consuming 
as the previous three. One can minimize this effort us- 
ing the exponential convergence method suggested and 
applied in Refs. [30. - .32.j . and/or the recently developed 
projected configuration interaction method |33l [3^ . In 
reverse, one can envision using some knowledge about 
the level density to extract the g.s. energy. This idea is 
not new (see, for example, Refs. [23l [35l - [37| ). However, 
we propose new algorithm that extracts the g.s. energy 
for a large model space by comparing the level density 
with that obtained in a reduced model space that can be 
exactly solved. 

The techniques described in this article are based on 
nuclear statistical spectroscopy 23]. We calculate the 
configuration spin and parity projected moments of the 
nuclear shell model Hamiltonian, which can be further 
used to obtain an accurate description of the nuclear level 
density up to about 15 MeV excitation energy. There- 
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fore, our methodology does not require any spin-cut-off 
parameter. One should mention that some of the more 
recent Monte Carlo approaches for level densities can also 
use direct spin projection techniques fl^]. 

The article is organized as follows. In Sec. II the 
fixed spin- and parity-dependent configuration moments 
method is revisited. The method allows one to trace 
such quantum numbers as parity and angular momen- 
tum explicitly. The extension of the algorithm to the 
proton-neutron formalism is discussed in Sec. III. Sec- 
tion IV is devoted to the results of the moments method, 
which are compared to exact shell model results and the 
results from Hartree-Fock-Bogoliubov plus combinatorial 
method. In Sec. IV we also present our new algorithm 
to extract the g.s. energy by comparing level density in 
related model spaces. Section V is devoted to conclusions 
and future prospects of the moments method. 



II. SPIN- AND PARITY-DEPENDENT 
CONFIGURATION MOMENTS METHOD 

In this work we closely follow the approach proposed in 
Refs. [H,!!!]. According to this approach one can calcu- 
late the nuclear level density p as a function of excitation 
energy E in the following way: 



(1) 



Here, a = {n, J,Tz,n} is a set of quantum numbers, 
where n is the number of particles (protons and neu- 
trons), J is total spin, Tj is isospin projection, and tt is 
parity, k represents a configuration of n particles dis- 
tributed over q spherical single-particle orbitals. Each 
configuration k is fixed by a set of occupation numbers 
K = {ki,K2,... where Kj is the number of parti- 

cles occupying the spherical single-particle level j. Each 
configuration has a certain number of particles, isospin 
projection, and parity. The sum in Eq. ([T]) spans all 
possible configurations corresponding to the given values 
of n^Tz, and tt. The dimension DaK equals the number 
of many-body states with given J that can be built for a 
given configuration k. The function Gan is a finite-range 
Gaussian defined as in Ref.|25]: 



Gqk(-E) = G{E + Egs — E, 
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(2) 
(3) 



where Ea,^ and aan are the fixed- J centroids and widths, 
which will be defined later. Egg is the ground state en- 
ergy, 77 is the cut-off parameter, and N is the normal- 
ization factor corresponding to the following condition: 
J-^ G(a:, a)dx = 1. In this work we treat as a free pa- 
rameter. From previous works (see for example (22| ) we 
know that the cut-off parameter is 77 3. We can slightly 
variate the value of 77 to achieve a better description of 
the nuclear level density, see Fig. 



The J-dependent moments method provides a good de- 
scription of the exact J-dependent shell-model level den- 
sity. Figure [T] presents the results for ^®Si in the sd-shell 
for different values of spin J and positive parity. Figures 
S] and [5] present results for ^^Fe and ^^Cr nuclei in the 
p/-shell. Similar results were obtained for the density of 
states using the general moments method (see examples 
in Refs. [22, 23, 38i]). A very important ingredient for 
our method is an accurate knowledge of the ground state 
energy Egg. It is also important to investigate the sensi- 
tivity of the results to the the cut-off parameter 77, and 
find optimal values for it. These issues will be discussed 
in more detail in Sec. IV. 

Let us define now the fixed- J centroids and widths 
from Eq. To calculate them for a two-body Hamil- 
tonian. 
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one has to calculate traces of the first and second power 
of this Hamiltonian, Tr[iJ] and Tr[i/^], for each configu- 
ration k: 



where 



(^)a«-Tr("'^^[^]/^-, 



(5) 
(6) 

(7) 
(8) 



Here the symbol of trace Tr'^"'"-' [• • •] means the sum of all 
diagonal matrix elements, as {i^, J\ ■ ■ ■ \v, J), over all 
many-body states {i^, J) belonging to the given configu- 
ration K and having a certain set of quantum numbers 
a, including spin J. Technically, it is more convenient 
to derive these traces in a basis with a fixed spin pro- 
jection Mz), Tr*-*'^^-'[- • •], rather than in the basis with 
fixed total spin J), Tr^"^-* [• • •]. J-traces can be eas- 
ily expressed through the Mz-traces, given the spherical 
symmetry of the Hamiltonian: 



Tr("')[---] = Tr' 



Tr 



(M,) 



(9) 



For simplicity, in Eq. (j9|) we omitted all quantum num- 
bers, except the projection Mz and the total spin J. 

Hereafter we use the label a to denote a set of quantum 
numbers that includes either the fixed Mz or the fixed J, 
keeping in mind that Eq. (j9|) can always connect them. 
In every important case we will point out which set of 
quantum numbers was used. Following the approach of 
Ref. t39i| (a similar method can be found in Ref. i40j ) , 
we obtain the following expressions for the traces from 
Eqs. and dH): 



(10) 
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FIG. 1: ^*Si, parity=+l. Comparison of nuclear level densities between exact shell model (solid line) and moments method 
(dashed line). Cut-oflF parameter rj = 2.8, interaction: USD, sd-shell. 
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where i,j,l, and q are single-particle states with cer- 
tain spin projections and possible occupation numbers 
equal to or 1. Notice that the single-particle orbitals 
we used to define the configurations in Eq. ([T|), can 
host all particles with all possible spin projections cor- 
responding to the orbital's spin. D^L = Tr^"''^ [a|ai] 
can be interpreted as a number of many-body states 
with fixed projection (if we consider M^-traces) and 
the single-particle state i occupied, which can be con- 
structed for the configuration k, Dai} = Tt^°"^'' [ala^jajai], 

and so on. These D- 



[alaj^a'ljaqajai], 



structures were called propagation functions in Ref. [33] . 
For completeness, we repeat here the recipe used to cal- 
culate them. One can show [soj that 



S<t <7l 



^ E 



(12) 



where all U are integers, configuration k' = 
{k[, ... , K^} can be derived from the original configu- 
ration K — K2, ... , Kq} by rcmoving t particles corre- 
sponding to the single-particle states ri, r2, • ■ ■ r^. Formal 
expression for the new k' configuration can be written as 
follows: 



E 



(13) 



where the sum includes only those values of i for which 
the corresponding single-particle state belongs to the 
single-particle level j. We also assume that all the oc- 
cupation numbers k'^ must be positive, which imposes 
certain restrictions on the possible values of the ampli- 
tudes ti. For every new configuration k' one can easily 
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FIG. 2: Speedup is defined as T\/Tn, where T„ is the cal- 
culation time, when n processors were used. These calcu- 
lations were performed on FRANKLIN supercomputer at 
the National Energy Research Scientific Computing Center 
(NERSC) [H. 



define new quantum numbers, a' = {n' M'^T!,tt'}, enter- 
ing Eq. (|12p . Examples are the new number of particles 
n' — n — t and the new spin projection: 



r-2 tsrur^, (14) 

where m,.. is the projection of the single-particle state 
Ti . The new isospin and parity tt' are defined similarly. 

III. THE MOMENTS METHOD ALGORITHM 
IN THE PROTON-NEUTRON FORMALISM 

Let us describe some technical features of the algo- 
rithm we developed for the nuclear level density calcu- 
lation. First of all, we treat protons and neutrons sepa- 
rately, that is, the basis of many-body wave functions are 
represented by a product of proton and neutron parts: 



|;.,iV4) = |i.p,MW).|^.„,Mi")), 



(15) 



where M^P^ + M^"' = M^. Thus, the wave functions 
(|15p have fixed isospin projection T^, but do not have 
a fixed isospin T. As we already mentioned, it is more 
convenient to use the basis of wave functions with fixed 
spin projection M^, rather than one with fixed spin J. 

One can gain essential advantages from such a proton- 
neutron separation of the basis. One of them is con- 
nected to the number of configurations that appear in 
the sum of Eq. ((T|). Naturally, the number of configu- 
rations with fixed is much greater than the number 
of configurations with fixed isospin. The large number 
of configurations allows the use of many-cores computers 
with greater efficiency. In other words, the calculation of 
the sum in Eq. ([T]) with a larger number of configurations 
can be more efficiently distributed on a larger number of 
processors. Figure [2] presents the speedup (calculation 



r(») 



speed gain) as a function of the number of used proces- 
sors. One can see that the case with the larger number of 
configurations, ^®Se, scales better than the case with the 
lower number of configurations, ^'^Ge. Up to 2000 cores 
the speedup is almost perfect (the dotted line presents 
an ideal speedup). At this point the calculation time is 
about of 1-2 minutes and further improvement is hardly 
achieved. 

Another significant advantage of the proton-neutron 
formalism is the new algorithm of calculation of the di- 
mensions Ua/t, -DqL, D^aK, and so on. Because of the 
proton-neutron separation one can calculate all proton 
and neutron dimensions separately. Later, the dimen- 
sions we are interested in can be easily constructed from 
the proton and neutron parts using the following convo- 
lution 
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(16) 



where, instead of the whole set of quantum numbers a, 
only the spin projection Mz was printed out. Kp and 
Kn are the proton and neutron parts of the configura- 
tion K. Equation pop can be easily applied to all types 

of dimensions, Da,\, we have shown in the formalism of 
Sec. II. The advantage comes from the fact that one 
can calculate and keep in memory all proton and neu- 



tron dimensions, D 



Mi' 



projections Mz and Mz 



and Z?^,(„,^ 



for all possible 



and for all possible configu- 
rations Kp and Kn- Afterwards, using Eqs. and p^ . 
one can very quickly calculate all the dimensions: -Dqk, 
-DqL, DaK, and so on, for all Mz and J. 

One more technical detail, which allows a significant 
speed up of the algorithm, is that using the proton neu- 
tron separation one can avoid multiple computations of 
the most time consuming structures, such as Dq^''''. Let 
us consider a case when all four single-particle states 
{ijql} are protons. One can then use an equation similar 
to Eq. (HH): 



D 
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D 
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(17) 



M'f'' + M'i"' = M, 



For all configurations k that have the same proton parts 



one will have to recalculate D 



Mi- 



for each neutron 

ilijqi] 



configuration. Alternatively, one can calculate D -^^Ip) 

only once, and store them in memory. That strategy, 
however, will require a large amount of storage. More 
efficiently, one can only store the contributions of the 
Dok'^^ structures to the width, Eq. pT|) . that is, one can 
only store the following structures. 
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where all single-particle states are protons. Thus, instead 
of using Eq. P7|) one can calculate the contribution to 
the width directly via the convolution, 



E 



(19) 

which is very similar to Eqs. ()16p and (|17p . As one 
can see the new approach avoids multiple calculations of 



D 



Storing the structures Eq. (|18p may signifi- 
cantly speed up the algorithm for large cases, such as 
^*Se in pf + 59/2 model space. The downside is that 
the calculation of the structures T, ,(p) ,T, does 
not always scale well on a large number of cores, since 
the number of these structures is much smaller than the 
total number of configurations. 



Element 


Space 


Total dim 


Elapsed time (sec) 


™Br 


pf + 99/2 




1.07 • 10* 


««Se 


pf + 99/2 


W'" 


1.03 • 10* 




pf + 99/2 




0.76 ■ 10* 


eoZn 


Pf 


10" 


37.4 


52Fe 


pf 


10"' 


13.6 


28Si 


sd 


10« 


0.7 



TABLE I: Elapsed times of nuclear level density calculations 
(for all J, positive parity) with the moments method code. 
The calculations were done on a 16 cores machine with 2.8 
GHz CPU frequency. 



Table U presents calculation times for different nuclei 
calculated in different shell-model spaces. The calcula- 
tions were done on a 16 cores machine with 2.8 GHz CPU 
frequency. One core ("master") distributed all the work 
between the other 15 cores ("slaves"). One can empha- 
size here that the listed times correspond to calculations 
of the nuclear level densities for all J and for positive par- 
ity. For the case of ^^Se the largest m-scheme dimension 
is about 10"'^^. For each J the m-scheme dimensions vary 
from 10^^ to 10^*, which makes direct diagonalization im- 
possible. Using the moments method and our algorithm 
we are able to calculate the shapes of nuclear densities 
for ^®Se in less than three hours on a 16 cores machine. 
If the number of processors reaches 1000 then one needs 
only a few minutes to complete the calculation. 



IV. RESULTS 

As a first example we consider the nuclear level den- 
sities of ^^Si in the sc?-shell model space, for which we 
use the USD interaction [4l|. Figure [T] presents the com- 
parison of the exact nuclear level densities of different 
spins (solid lines) with those obtained with the moments 




1,8 
1.6 
1,4 
, S 1,2 

% > 
g 0,8 

0,6 

0,4 

0,2 



2,5 3 3.5 

cut-off parameter r[ 

FIG. 3: Error factor /err as function of the cut-off parameter 



methods (dashed lines). Equations ([TJ and ^ require 
the knowledge of the ground state energy Egg and the 
cut-off parameter rj. While the ground state energy of 
^*Si can be calculated in this case using the standard 
sheU model. Egg = -135.94 MeV, for the value of the 
cut off parameter r] we only have some general idea that 
it should be around 3 [l^ HBl • In Fig- HI by choosing 
T] — 2.8, the moments method reproduces quite well the 
exact shell model densities. To get a better description of 
the moments method level densities we can adjust the rj 
parameter to optimally reproduce the exact shell-model 
densities. The cut-off parameter plays a similar role as 
that of the width in a Gaussian distribution. Indeed, if 
we increase the cut-off parameter then the density be- 
comes wider and lower, while decreasing it leads to a 
narrowing of the density. Fig. |3] helps to determine the 
optimal value of the cut-off parameter ry. In this figure 
the vertical axis presents an error factor ferr, which is a 
measure of the deviations of the calculated density pmm 
(using the moments method) from the exact shell-model 
level density, psm- One possible way to construct this 
error factor is the following [Z^, ■ 



ferr — exp 



i=l 



- 1, (20) 



where the sum over i spans an energy region, for which 
one wants to compare the level densities. The moments 
method is known to be statistically valid when the fluc- 
tuations can be described by a Gaussian orthogonal En- 
semble. Strictly imposing this condition may not be very 
practical. Therefore we consider it valid in the regime 
where the density is at least 5-10 levels per MeV. This 
condition can be used to establish a starting energy for 
the sum in Eq. pO|. The sum should also be upper 
limited to excitation energies for which the 2?iuj states 
are not contributing significantly. For most of the model 
spaces considered here this upper limit is about 10-15 
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FIG. 4: ^^Fe, parity=+l. Comparison of nuclear level densities between exact shell model (solid line) and moments method 
(dashed line). Cut-off parameter rj = 2.6, interaction: GXPFIA, p/-shell. 



MeV. Figure |3] presents the dependence of the ferr on 
the cut-ofT parameter rj. It suggests that optimal values 
for 77 are in the interval 2.5-3.0, which supports our initial 
guess. It also indicates that there is relatively small sensi- 
tivity to this parameter in the indicated interval. There- 
fore, for the pf and pf + gg/2 spaces we chose rj = 2.6, 
the value for which the moments method level densities 
reproduce quite well the exact shell-model level densi- 
ties. The minimum value of ferr — 0.1 — 0.2 offers an 
estimated average accuracy of the moments method for 
the model spaces and the nuclei shown in the inset. A 
study of an optimal rj parameters for a larger class of 
nuclei and model spaces will be published elsewhere [i^ . 
One should also mention that the exact spin- and parity- 
dependent shell-model densities were calculated with the 
NuShellX code H. 

Next we present a couple of examples for the »/ shell, 
for which we used the GXPFIA interaction [4^[4^. Fig- 
ures S] and [S] present the results for ^^Fe (J — 0,1) and 
for ^^Cr (J = 0, 1, 2, 3) in the p f sh ell , for which we have 
used the GXPFIA interaction [45l. l46l| . The correspond- 
ing ground state energies are known and the cut-off pa- 
rameter was chosen to be 77 = 2.6. One can only compare 
the lowest parts of the level densities (up to 200 levels). 
For higher excitation energies it already becomes too dif- 
ficult to calculate the exact shell-model densities, because 
of the large number of states needed. As one can see, the 
moments method densities are in a very good agreement 
with the exact shell-model densities. Figures |4] and [5] 
also includes results obtained by Goriely, et al. using the 
Hartree-Fock-Bogoliubov (HFB) sing le particle energies 
and the combinatorial method [3, |43] . 

We mentioned in the introduction that one can envi- 
sion using information from level densities to extract with 
good approximation values for the g.s. energies. Using 
our algorithm and the moments method one can easily 
calculate the nuclear level density for any nucleus that 
can be described in the pf + 59 /2 model space. The in- 
teraction we used for this model space was built starting 



with the GXPFIA interaction for the pf model space, 
to which G-matrix elements that describe the interac- 
tion between the pf orbits and 39/2 orbit were added. 
The single particle energy for the (79/2 orbit was fixed at 
-0.637 MeV. The calculation time for the worst case is 
more than reasonable: it takes about three hours for 16 
processors and only a few minutes for 1000 processors. 
Figures [5] and [7] summarize the results obtained for ^^Se 
and ^''Ge, nuclei that are believed to be "waiting-points" 
in the rp-process path (48l - [50| . We have only presented 
the densities for J = 0, 2, and positive parity. 

It is important to notice that in the pf model space 
the shell-model calculations of the g.s. energies can be 
done. For the pf shell we have the following ground 
state energies: Egs{pf) = -304.25 MeV for ^^Ge and 
Egsipf) = -353.1 MeV for ^SSe. Using these ground 
state energies and the cut-off parameter 77 = 2.6, we are 
able to calculated the nuclear level densities according to 
Eqs. (dJ and The solid lines on Figs. [S] and [7] present 
the densities in the pf shell. To calculate the same level 
densities in the pf + gQ/2 model space we have to adjust 
the ground state energies and the cut-off parameter for 
this space. For the cut-off parameter we use the same 
value, T] — 2.6. The dotted lines show the nuclear level 
densities if we keep the ground state energies as they 
were in the pf shell. It is natural to expect only small 
differences between the level densities calculated in those 
two model spaces at low excitation energy. The differ- 
ence must be compensated by the fact that the ground 
state energy for the larger model space, that is pf + (79/2, 
must be lower compared to the ground state energy for 
the smaller model space, that is, pf. By decreasing the 
ground state energies for the pf + gQ/2 model space, one 
gets the dashed lines on Figs. [S]and[71 The dash-dotted 
lines on Fig. [S] correspond to an ground state energy 
Egs = -305.8 MeV of ^^Ge, which was obtained by a 
truncated shell model calculation in which up to six par- 
ticles were excited from the /7/2 orbits and/or into the 
(79/2 orbit. The 7n-scheme dimension in this calculation. 





FIG. 5: ^^Cr, parity=+l. Comparison of nuclear level densities between exact shell model (solid line) and moments method 
(dashed line). Cut-off parameter rj = 2.6, interaction: GXPFIA, p/-shell. 




FIG. 6: Ge nuclear level densities for J = 0, 2 and positive parity. Solid line presents the calculation in the pf shell with 
GXPFIA interaction. For this calculation wc know the ground state energy Egsipf) = —304.25 MeV. Other three lines present 
calculations in the large model space, when level 39/2 is added. The ground state energy for these cases Egs{pf + (79/2) = 
Egs{pf) — AE, where AE is the energy shift. The cut-off parameter is 77 = 2.6. 




Excitation energy (MeV) Excitation energy (MeV) 



FIG. 7: Se nuclear level densities for J = 0, 2 and positive parity. Solid line presents the calculation in the pf shell with 
GXPFIA interaction. For this calculation we know the ground state energy Egs{pf) — —353.1 MeV. Other two lines present 
calculations in the large model space, when level gg/2 is added. The ground state energy for these cases Egs{pf + gg/2) = 
Egsijpf) — Ai?, where Ai? is the energy shift. The cut-off parameter is 77 = 2.6. 



13.5 X 10^, is at the upper limit of the state-of-the-art 
shell model calculation. As one can see, this value does 
not describe satisfactorily the low excitation energy level 
densities. In order to make the low- lying part of the two 
densities very close, one has to adjust the ground state 
energies for the pf -\- 179/2 model space to the following 
values: 

Egs{pf + 59/2) = -306.7 MeV for Qe, (21) 

Egs{pf + 59/2) = -356.5 MeV for Se. (22) 

The " low- lying part of the density" should be chosen such 
that the excitations to the ji orbit do not give a signifi- 
cant contribution. For these cases we use the interval 3-6 
MeV in excitation energy. We conclude that the g.s. en- 
ergy adjustment of Eqs. (pij) and (|^ can be treated as a 
method of estimating the ground state energies in larger 
spaces. Therefore, one can formulate now the following 
recipe: to get the ground state energy for a nucleus in a 
large model space, in which the direct shell model calcula- 
tion is presently impossible, one can calculate the nuclear 
level densities in the large model space and in an asso- 
ciated smaller model space, for which the ground state 
energy calculation is possible. Then, the ground state 
energy for the larger model space can be estimated by de- 
manding that the level densities in the two model spaces 
at low excitation energy be the same or very close. Cer- 
tainly, one should not arbitrarily select the larger model 
space and the associated solvable model space. What we 
proved here seems to be valid when adding one more sin- 
gle particle level to a solvable model space, such that the 
entire shell structure is not significantly distorted. 

V. CONCLUSIONS AND OUTLOOK 

In conclusion, we developed a very efficient algorithm 
to calculate the shell-model spin- and parity-dependent 



configurations centroids and widths, which can be used 
to calculate nuclear level densities. The new algorithm 
takes advantage of the separation of the model space in 
neutron and proton subspaces. This separation provides 
two important advantages: (i) the exponentially explod- 
ing dimensions and propagators can be calculated more 
efhciently in proton and neutron subspaces, and the full 
results can be recovered via simple convolutions; (ii) the 
number of configurations is significantly increased in the 
proton-neutron formalism, which very much improves the 
scalability of the algorithm on massively parallel comput- 
ers. Our tests indicate almost perfect scaling for up to 
4000 cores, and we are convinced that it can scale well 
up to tens of thousand of cores. The new algorithm is 
so fast that the bottleneck of the calculation is now that 
of the ground state energy. That is why we cannot test 
our algorithm for cases that take more than 1 minute on 
4000 cores. 

Therefore, we investigated the possibility of using the 
calculated shapes of the nuclear level densities to extract 
the g.s. energy. We showed that by slightly increment- 
ing the model space, and imposing the condition that 
the level density does not change at low expectation en- 
ergy, one can reliably predict the g.s. energy, and further 
the full level density. This new method of extracting the 
shell model ground state energy for model spaces whose 
dimensions are unmanageable to direct diagonalization 
opens new opportunities for calculating shell model nu- 
clear level densities of heavier nuclei of interest for nuclear 
astrophysics, and nuclear energy and medical physics ap- 
plications. In particular, one can envision using effec- 
tive interactions extracted from ab-initio theories, such 
as the G-matrix with core polarization, with some ad- 
justable monopole corrections that can be tuned to de- 
scribe the effect of the correlations to nuclear level den- 
sities of heavy nuclei. This class of effective interactions 
is much larger than the class of pairing plus multipole 
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interactions that Monte Carlo methods jl5| can use. 

The present method could be also used in more than 
one major harmonic oscillator shell for medium-mass nu- 
clei to describe level densities of both parities. The 
center-of-mass spurious states could be eliminated in 
these cases using a method proposed in Rcf. [27|. This 
method requires an extension of the present algorithm 
that will enforce restrictions on the classes of configu- 
rations included in the widths formula, similar to the 
one proposed in Ref. [1^. Work in this direction is in 
progress. 

Our method seems to exhibit some sensitivity to the 
cut-off parameter -q. In the cases we studied a value of 
about 2.8 seems to provide very good results, but further 
investigations of the optimal values of this parameter are 
necessary. In addition, one should consider going beyond 
the two-moments approach for the configuration distri- 



butions. These higher moments were used in the past for 
the density of states. The J-dependent higher moments 
are more difficult to calculate, but given the computa- 
tional advances we made with the first two moments, 
one will envision an efficient algorithm to calculate the 
higher moments in the near future. 
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